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We compute the thermal axino production rate in supersymmetric QCD to lead- 
ing order in the gauge coupling. Using hard thermal loop resummation and the 
Braaten— Yuan prescription, we obtain a finite result in a gauge-invariant way, 
which takes into account Debye screening in the hot quark— gluon-squark—gluino 
plasma. The relic axino density from thermal reactions in the early Universe is 
evaluated assuming the axino is the lightest supersymmetric particle and stable 
due to R-parity conservation. From the comparison with the WMAP results, we 
find that axinos could provide the dominant part of cold dark matter, for example, 
for an axino mass of 100 keV and a reheating temperature of 10® GeV. 

1. Introduction 

In supersymmetric extensions of the standard model in which the strong 
CP problem is solved by the Peccei-Quinn (PQ) mechanism, the axino 5, 
arises naturally as the fermionic superpartner of the axion.^ Thus, the axino 
is electrically and color neutral and its interactions with the MSSM particles 
are suppressed by the PQ scale fa/N 5 x 10^ GeV. If the axino is the 
lightest supersymmetric particle (LSP) and if i?-parity is conserved, axinos 
are stable and could provide the dominant part of cold dark matter. 

In this talk we present the computation of the thermal axino production 
rate at high temperatures and discuss its cosmological implications. With 
inflation governing the earliest moments of the Universe, any initial pop- 
ulation of axinos was diluted away and the thermal production of axinos 
set in at the reheating temperature Tji. We restrict our investigation to 
fa/N > Tr > 10'' GeV, where the U(1)pq symmetry is broken and axino 
production from decays of particles out of equilibrium is negligible.^ The 
results presented are extracted from Ref. 3 where more details can be found. 



•Talk given by F.D.S. at SEWM 2004, Helsinki, Finland, June 16 - 19, 2004. 
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2. Thermal Axino Production in Supersymmetric QCD 

We concentrate on the axino-gluino-gluon interactions given by the 
diniension-5 interaction term 

with the strong gauge couphng g, the gluon field strength tensor GJ^j,, the 
gluino g, and A'' being the number of quarks with PQ charge. The resulting 
2 — > 2 scattering processes in supersymmetric QCD and the corresponding 
squared matrix elements are listed in Table 1, where s = {Pi + ^2)^ and 
t = (Pi — Ps)^ with Pi, P2, P3, and P referring to the particles in the given 
order. Working in the limit, T ^ rrii, the masses of all particles involved 
have been neglected. Sums over initial and final spins have been performed. 
For quarks and squarks, the contribution of a single chirality is given. 

The processes B, F, G, and H lead to a logarithmic coUinear singu- 
larity due to t-channcl and ?i-channel exchange of soft (massless) glu- 
ons. Here screening effects in the hot quark-gluon-squark-gluino plasma 
(QGSGP) have to be taken into account. With hard thermal loop (HTL) 
rcsumination,'* this can be done systematically in a gauge-invariant way. 
Following the Braatcn Yuan prescription,^ we introduce a momentum scale 
fccut such that gT <C fccut <C T in the weak coupling limit, g 1. This sep- 
arates soft gluons with momentum transfer of order gT from hard gluons 
with momentum transfer of order T. In the region of soft momentum trans- 



Table 1. Squared matrix elements for axino (a) produc- 
tion in two-body processes involving quarks (qi), squarks 
((ji), gluons (g"), and gluinos (g") in the high-temperature 
limit, T > mi, with the SU(3) color matrices /"'^ and T?. 



process? \Mi\y j^^^^:ifj^j^ 
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9" 


+ 9'' -9" 


+ a 


B 


9" 


+ 9''-* 9'= 


+ a 


C 


qi 


+ 9''^qj 


+ a 


D 


ga 


+ 9i — qj 


+ d 


E 


qi 


+ qj-^9'' 


+ a 


F 


9" 


+ g''^g'' 


+ a 


G 


qi 


+ g''^qj 


+ a 


H 


qi 


+ 9''^qj 


+ a 


I 


qi 


^qj^r 


+ a 


J 


qi 


+ qj-*9'' 


+ a 



4{s + 2t + 2i^)\r'"'\'^ 
-4{t + 2s-f-2^)|/°''=|2 

-2t|r;j2 
-2t|r;j2 



st(s + t) 'J I 

-2(1 +2s + 24)|T^|2 

-4(t + 4)|r;j2 

2(1 +2t + 2^)|T?.|2 
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fer, k < fccut, we use the HTL-resummed gluon propagator, which takes into 
account Debye screening in QGSGP in terms of the supersymmetric ther- 
mal gluon mass = gT^/ [Nc + nf)/6, where iVc = 3 and n/ = 6 are 
respectively the number of colors and color triplet and anti-triplet chiral 
multiplets. Computing the imaginary part of the thermal axino sclf-cncrgy 
with the ultraviolet (UV) momentum cutoff fccut) the soft contribution to 
the production rate of axinos with energies E "^T is obtained 



= fF{E) 



3g\N^ - l)mlT 



In 



1.379 



(2) 



where /f(-E) = [exp(£^/T + In the region of hard momentum trans- 

fer, fc > fccut, the bare gluon propagator can be used. From the squared 
matrix elements of the 2^2 processes (cf. Table 1) weighted with ap- 
propriate multiplicities, statistical factors, and phase space densities, the 
corresponding hard contribution can then be obtained conveniently 



hard 



32 7r4 {fjNf 



In ' 



1287r4 



I'CUt 



(3) 



with fccut as the infrared (IR) cutoff for gluon momentum transfers. The ex- 
pression represented by the ellipses can be found in Ref. 3. It is independent 
of fccut- Thus, by summing the soft and hard contributions, the artificial 
fccut dependence cancels and the finite leading order rate is obtained 



(4) 



hard 



<Pp d?P soft 

In Fig. la, we show the normalized thermal axino production rate 
l/V-a dTa/diE/T) as a function of E/T for temperatures of T = 10*^ GeV, 
lO^GeV, lO^GeV, and lO^GeV. For E < T, we expect modifica- 
tions as our computation is restricted to the thermal production of hard 
{E ^ T) axinos.^ Here and below the 1-loop running of the strong cou- 
pling in the MSSM is taken into account by replacing g with g{T) = 
[g~^{Mz) + 31n{T/Mz)/i87r'^)]^^/^, where the value of the strong coupling 
at the Z-boson mass Mz, g^{Mz)/{4:TT) = 0.118, is used as input. For 
E/T < 2, the rate turns negative except for very high temperatures. With 
g{T= 10^ GcV) = 0.986 and .g(T= 10^ GcV) = 0.880, this unphysical be- 
havior follows from the extrapolation of our result from 5 <C 1 to <? « 1. 
As higher-order corrections might become sizable for T ^ 10^ GeV, new 
techniques are needed to compute the thermal production rate reliably also 
for g ^ I. At present, however, our result derived in a gauge-invariant way 
supersedes previous gauge-dependent and cutoff-dependent estimates.^ 
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Normalized Thermal Axino Production Rate Axino Density from Thermal Production 
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Figure 1. (a) Tlic normalized thermal axino produetion rate l/TadTa/d{E/T) as a 
function of E /T for T = 10^ GeV (dotted line), lO'^ GeV (dash-dotted hne), 10* GeV 
(daslied hne), and lO" GeV (solid line). The curves result from (4) derived for E ^ T. 
(b) The axino density parameter Clah as a function of Tr for fa/N = lO" GeV and 
ma = 1 keV (dotted line), 100 keV (dashed line), and 10 MeV (solid line). The grey 
band indicates the WMAP result on the relic density of cold dark matter (2cr error). 



3. Stable Axinos as Dark Matter in the Universe 

By comparing the axino interaction rate for /a/iV = 10^^ GeV with the 
Hubble parameter H for the early radiation-dominated epoch, an axino de- 
coupling temperature of Tp w 10^ GeV is estimated.^ For Tf{ < Tjj, axinos 
have not been in thermal equilibrium after inflation. Here the evolution 
of the axino number density na with cosmic time t is described by the 
Boltzmann equation with a collision term Cg accounting for both the axino 
production and disappearance processes in the primordial plasma, 

^ + 3Fna = Ca. (5) 

at 

The disappearance processes can be neglected for Tr sufficiently below Tq. 
Then, by integrating the production rate (4), we obtain the collision term 



^ {N^ - 1) 3C(3)g6T6 
"~ (/„/7V)2 40967r7 



(1 380 T"^ \ 
- — 2 — (Arc + n/) + 0.4336 n/ 
nig ) 



(6) 



where C(3) « 1.2021. Assuming conservation of entropy per comoving 
volume, the Boltzmann equation can be solved analytically. The resulting 
present (to) axino density parameter fJa^^ = iTiana{to)h'^ / pc with pc/h? = 
3.6 X 10~^GeV depends on the axino mass rria, the PQ scale fa/N, and 
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the reheating temperature Tr in the following way 

2 



where g = ^(Tr). In Fig. lb, this result is illustrated as a function of Tr 
for fa/N = 10" GeV and ma = 1 keV, 100 keV, and 10 MeV. For Tr 
above To, axinos are in thermal equilibrium so that Og/i^ is independent 
of Tji as shown for m,a = 1 keV by the horizontal line. There will be a 
smooth transition instead of a kink once the axino disappearance processes 
are taken into account. The grey band indicates the WMAP result'' on the 
cold dark matter density (2cr error) fig^^/i^ = 0.113lg:gJt 



4. Conclusion 

For ma = 100 keV and Tr 10^ GeV, the relic axino density (obtained 
with fa/N = 10" GeV) agrees with the WMAP result on the cold dark 
matter density. Although relatively light for being cold dark matter, axi- 
nos with ma = 100 kcV could still explain large-scale-structure formation, 
the corresponding power spectrum, and the early reionization observed by 
WMAP. Thus, as far as abundance and structure formation are concerned, 
axinos from thermal production in the early Universe arc indeed a viable 
solution of the cold dark matter problem.* This conclusion has already been 
drawn in Ref. 2 based on a gauge-dependent and cutoff-dependent estimate 
of the relic axino abundance from thermal production. Since our result for 
flah^ is smaller by a factor of ten, the upper limit on Tr for a given value 
of ma is relaxed by one order of magnitude. This can be important for 
models of inflation and for the understanding of the baryon asymmetry in 
the Universe. Already Tr « 10^ GeV is relatively small and excludes, for 
example, thermal leptogenesis. 
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